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ABSTRACT 

Radiative transfer studies of Type la supernovae (SNe la) hold the promise of constraining 
both the time-dependent density profile of the SN ejecta and its stratification by element 
abundance which, in turn, may discriminate between different explosion mechanisms and 
progenitor classes. Here we present a detailed analysis of Hubble Space Telescope ultraviolet 
(UV) and ground-based optical spectra and light curves of the SN la SN 2010jn (PTFlOygu). 
SN 2010jn was discovered by the Palomar Transient Factory (PTF) 15 days before maximum 
light, allowing us to secure a time-series of four UV spectra at epochs from — 1 1 to +5 days 
relative to S-band maximum. The photospheric UV spectra are excellent diagnostics of the 
iron-group abundances in the outer layers of the ejecta, particularly those at very early times. 
Using the method of 'Abundance Tomography' we have derived iron-group abundances in 
SN2010jn with a precision better than in any previously studied SNIa. Optimum fits to the 
data can be obtained if burned material is present even at high velocities, including significant 
mass fractions of iron-group elements. This is consistent with the slow decline rate (or high 
'stretch') of the light curve of SN 2010jn, and consistent with the results of delayed-detonation 
models. Early-phase UV spectra and detailed time-dependent series of further SNe la offer a 
promising probe of the nature of the SN la mechanism. 

Key words: supernovae: general - supernovae: individual (SN2010jn) - techniques: spec- 
troscopic - radiative transfer 



1 INTRODUCTION 

Supernovae (SNe) play an important role in many areas of modern 
astrophysics. In particular, Type la SNe (SNe la) produce most of 
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the iron-group elements in the cosmos l llwamoto etahll 19991) . and 
can be used as 'standardi zable candles' to probe the expansion his- 
tory of the unive rse (e.g. iRiess et alJll998l;|Perlmutter et aljl999t 



iRiess et al.l2007l : lKessler et alj2009l : lsullivan et alj201lh . Because 
of their importance, extensive efforts are under way to build a 
comprehensive observational and theoretical picture of SNela (e. 
Hillebrandt & NiemeveJl2000l: iMazzali et al]|2007l : lHowell|[201 



Ropk e et al EoiTh . A primary aim is to understand the nature and 



diversity of the explosions and thus to place on a physical ground 
the empirical calibration procedures that are used to deduce the lu- 
minosity of SNe la f rom their light-curve shape (e.g. IPhillipsI 19931 ; 
|Phillipsetal.1ll999l) . Another aim is to determine the progenitor 
systems of SNe la. 

One way to improve our understanding of SN la physics is to 
analyse and interpret detailed observations through radiation trans- 
port models. A particularly useful technique is that of 'Abundance 
Tomography' dStehle et alj I2005T) . where a time-series of SNIa 
spectra is modelled and information obtained about the density pro- 
file of the SN ejecta and its stratification in terms of element abun- 
dances. This makes it possible to describe the mode of explosion 
and possibly to discriminate among different progenitor scenarios. 
Many detailed spectral time series for SNela have become avail- 
able in the last decade, but most of them are restricted to optical 
wavelengths. The ultraviolet (UV) spectrum of SNela is shaped 
by iron-group elements, which dom i nate line-block ing and fluores- 
cence effects jPauldrach et alJl996l : lMaz zali 2000). Thus, series of 
UV spectra are invaluable diagnostics for iron-group abundances 
and ejecta densities in different zones of the SN. These in turn are 
indicators for basic e xplosion properties s uch as burning efficiency 
or explosion energy dlwamoto et alj[l999l) . Extensive work on UV 
spectra of SNe la has been carried out in the last few years (e.g. 
Ellis et al.ll2008l; iFolev et alj|2008l: Iskuer et alJ[20o3 : ICooke et al.1 



201 ll : IWalker et alj|2012t iMaguire et alj|2012h . however the focus 



has usually been to obtain a single spectrum per SN. Both ob- 
servations and models indicate that even SNel a with similar op- 
tical spectra can be very different in the UV dLentzetalJl2000l : 
lEllis et al.l20b"&ISauer et alj|2008l : lwarker et alj2012h . 

Here we present observations and models of the SNIa 2010jn 
(PTFlOygu). The SN was discovered by the Palomar Transient Fac- 
tory (PTF) only a few days after explosion, which allowed a de- 
tailed spectral time series to be obtained, including UV data from 
the Hubble Space Telescope. The combined optical-UV spectra 
make it possible to conduct an analysis of the SN in unprecedented 
detail. Abundance Tomography was performed based on density 
profiles resulting from single-degenerate Chandrasekhar-mass ex- 
plosio n models. These models (cf. e.g. iHillebrandt & NiemeveJ 
120001) assume that SNela are the explosion of an accreting 
carbon-oxygen (CO) w hite dwarf (WD). We used both the fast - 
deflagration model W7 dNomoto et al.lll984llwamoto et ai]|l999l). 
which has been shown to match average SNe la Jstehle et al] 
l2005l : lTanaka et alj201 l]), and a more ener getic delayed-detonation 
model (WS15DD3. Ilwamoto et al.1 Il999h . The latter model as- 
sumes that the combustion flame, which initially propagates sub- 
sonic ally (deflagration ), becomes supersonic (detonation) at some 
point dKhokhlovll l_99l|). Based on the quality of the fits, we suggest 
that the latter model is a more realistic description of SN 2010jn. 

The paper starts with a report on the observations and a pre- 
sentation of the observed spectra, which were obtained from —13 d 
to +5 d (time in the SN rest frame) with respect to maximum light 
in the rest-frame B-band (Section[2]l. Afterwards, we discuss our 
models. We lay out the objectives, methods and assumptions (Sec- 
tions [3T|-[3]4](. Synthetic spectra are presented and the inferred 



abundance profiles are discussed (Sections l3, 5113,61 ). Finally, results 
are summarised and conclusions are drawn (Section[4]l. 



2 OBSERVATIONS, DATA REDUCTION AND 
PHOTOMETRIC PROPERTIES 

SN2010jn (PTFlOygu) was discovere d on 2010 October 12 UT 
by th e Palomar Transient Factory (PTF: IL"aw et al . 2009; R au et al.1 
120091) using the Palomar 48-in telescope (P4 8), through the ci ti- 
zen science project 'Galaxy Zoo: Supernova' dSmithet allboillfl 
The SN was found at a magnitude of r ~ 19.2 in the Sbc galaxy 
NGC2929, at a heliocentric redshift of z = 0.02505 (CMB-frame 
redshift of 0.02602, yielding a distance modulus fj, — 35.2 majj3), 
and a spatial position of RA 09:37:30.3, Dec. +23:09:33 (J2000). 
Only a marginal detection was present in data from the previous 
night, and no detection in data taken on 2010 October 08. Because 
of the relatively low redshift and apparent early discovery, a classi- 
fication spectrum was triggered using the Gemi ni-N telescope and 
the Gemini Multi-Object Spectrograph (GMOS: lHook 



lescope and 
et al.ll2004 



on 2010 October 13 (program ID GN-2010B-Q-13). This revealed 
an early SNIa at about 15 days before maximum light. 

Based on this classification, four epochs of UV observations 
of this SN were triggered with the Hubble Space Telescope (HST) 
using the Space Telescope Imaging Spectrograph (STIS) as part of 
the cycle 18 program 12298: 'Towards a Physical Understanding of 
the Diversity of Type la Supernovae' (PI: Ellis). A spectral monitor- 
ing campaign with ground-based telescopes was also commenced. 
Multi-colo ur light curves we re obtained with the Liverpool Tele- 
scope (LT; ISteele et alj[2004h using RATCam and with the Faulkes 
Telescope North (FTN). 

An observing log for the HST data can be found in TableQ] We 
use two observing modes, both with a 0.2" slit: the G430L/CCD on 
all four epochs (giving coverage from ~2900 A out to ~5700 A), 
and the G230LB/CCD (with nominal coverage from ~2000 A to 
~3000 A). Unfortunately, the more sensitive G230L/MAMA mode 
was not available as the MAMA was offline. Since the intrinsic UV 
flux of SNela below ~2700 A is low, the G230LB/CCD observa- 
tions only provide useful data redwards of this wavelength. 

We also make use of five ground-based optical spectra, taken 
with Gemini-N with GMOS and the William Herschel Telescope 
(WHT) with the Intermediate dispersion Spectrograph and Imag- 
ing System (ISIS). With GMOS, we used the B600 grating with a 
central wavelength of 450nm, and the R400 grating with a central 
wavelength of 750nm. With ISIS, we used the R158R (red arm) and 
R300B (blue arm) gratings, together with the 5300 dichroic. Again, 
details can be found in Table Q] At the earliest epoch we anal- 
yse (— 12.9 d), only an optical spectrum is available. At all other 
epochs, we use both HST and ground-based spectra, taken within 
12 hours of each other. 



2.1 Data reduction 

2.1.1 Spectra 

The tr eatment of the HST spectra is described in iMaguire et al.1 
Briefly, the spectra were downloaded from the HST archive 



1 PTFlOygu discoverers: Peter Woolliams, tracey, Graham Dungworth, lp- 
spieler, John P Langridge, Elisabeth Baeten, Tomas Raudys, adam elbourne, 
Robert Gagliano 

2 A Hubble constant Hq = 72kms _1 Mpc~ 1 is used throughout. 
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Table 1. Log of the spectroscopic observations of SN2010jn/PTF10ygu. 



Calendar date of 
observation (UT) 


Date of 
obs. (MJD) 


Phase 
(days)" 


Telescope 


Instrument 
configuration 


Exposure 
time (s) 


Oct 13 


55482.6 


-12.9 


Gemini-N 


GMOS/B600/R400 


450 


Oct 16 


JJ^OJ .\j 


— 10 5 


HST 




2175 


Oct 16 


55485.6 


-10.0 


Gemini-N 


GMOS/B600/R400 


450 


Oct 20 


55489.8 


-5.9 


HST 


STIS/G430L/CCD 


2175 


Oct 20 


55489.9 


-5.8 


HST 


STIS/G230LB/CCD 


4920 


Oct 21 


55490.6 


-5.1 


Gemini-N 


GMOS/B600/R400 


450 


Oct 26 


55495.4 


-0.4 


HST 


STIS/G430L/CCD 


2175 


Oct 26 


55495.5 


-0.3 


HST 


STIS/G230LB/CCD 


7665 


Oct 27 


55496.6 


+0.8 


Gemini-N 


GMOS/B600/R400 


450 


Oct 31 


55500.7 


+4.8 


HST 


STIS/G430L/CCD 


2175 


Nov 01 


55501.2 


+5.3 


WHT 


ISIS/R300B/R158R 


900 



days in the SN rest frame relative to maximum light in the rest-frame B-band 



using the on-the-fly reprocessing (OTFR) pipeline, giving fully cal- 
ibrated and extracted 1-D spectra, where the reduction and extrac- 
tion is optimised for point sources. The OTFR pipeline uses the 
latest calibration files and data parameters to perform initial 2-D 
image reduction such as image trimming, bias and dark current 
subtraction, cosmic-ray rejection, and flat-fielding. It then performs 
wavelength and flux calibrations. We applied further cosmic-ray re- 
moval and hot-pixel masking by hand. The ground-based spectra 
were reduced using standard IRAF procedures (see acknowledge- 
ments), including bias subtraction, flat-fielding, wavelength cali- 
bration, flux calibration, and telluric feature removal. The te ch- 
niques used were similar to those described in lEllis etalJi2008h . 

For all epochs except the first one we combined the various 
spectra available (i.e., one or two HST spectra and one ground- 
based spectrum). We used HST data below ~ 5500 A, and ground- 
based data above. In each case, we matched the spectra in the over- 
lap region using the HST spectrum as reference. The resulting four 
optical /UV spectra analysed here sample the SN in regular time 
intervals of ~5 day s. They are publically a vailable through the 
WISeREP archivtFI jYaron & Gal-Yamll2012l) . We have corrected 
all final 1-D spectra for Milky Way extinction using Ry — 3.1, 
a colour exces s of E (B — V)mw = 0.03 from the dust maps of 
Schlegel et alj dl998l). and t he M ilky Way dust extinction law of 
Cardelli. Clayton & Mathisl dl989i) . 

2.1.2 Light curves 

Light curves were produced from the P48 R and LT/FTN gri pho- 
to metric data using diff erence imaging; full details can be found 
in iMaguire et all j2012h . In all cases, a reference image with no 
SN light is produced, and subtracted from each image in which the 
SN is present by matching the point-spread functions (PSFs) of the 
images. In the case of the P48 data, this reference image was con- 
structed using data from before the SN exploded; for the LT/FTN 
data the reference image was taken 18 months after the SN was 
first detected, by which time the SN flux was < 1% of its peak 
value. We measure the SN photometry using a PSF fitting method. 
In each image frame (prior to subtraction), the PSF is determined 
from nearby field stars. This average PSF is then fit at the position 
of the SN in the difference image. Each pixel is weighed according 
to Poisson statistics, yielding a SN flux and flux error. 

The light curves are flux-calibrated to the Sloan Digital Sky 

3 |http : / /www . weizmann . ac . il/astrophysics/wiserep/| 



Survey (SDSS) D ata Release 8 (DR8 ) photometric system, close 
to the AB system dOke&GunrJl 19831) . using stars in the field-of- 
view of the SN. We generally follow the procedures outlined in 
lOfek et al.1 j2012h . In the LT/FTN data (g, r and i filters), the colour 
terms to the SDSS system are very small. The P48 colour term is 
larger, and we include colour, airmass, and colour-airmass terms. 
The r.m.s. of the colour term fits is ~0.02 mag (with a colour term 
in r — i of ~0.22 mag), although we present magnitudes in the natu- 
ral P48 system and do not apply this colour term to our photometry. 
Finally, we fit the same PSF used for the SN flux measurement to 
the calibrating field stars. We then use these data to apply aper- 
ture corrections to the SN photometry on an image-by-image basis. 
These corrections are small (<3%) and account for the difference 
between the fixed aperture in which the calibration is performed, 
and the PSF-fit in which the SN photometry is measured. In Fig.Q] 
our observed light curve data are presented. 

2.2 Photometric properties 

Estimating the peak luminosity of SNela in a given rest-frame 
bandpass requires an interpolation between observed data points 
at the time of maximum light, follow ed by a fc-correction (e.g. 
iNugent et al]|2002l ; Irlsiao et all l2007t) back to the standard rest- 
frame filter of interest using a time series spectral energy distri- 
bution (SED). We fit the obser ved photometry for SN2010jn us- 
ing the SiFTO light curve fitter dConlev et alj|2008l) developed for 
SNIa cosmology studies. 

SiFTO works in flux space, manipulating a SED template and 
synthesizing an observer-frame light curve from a given spectral 
time-series in a set of filters at a given redshift. These template 
light-curves are then fit to the data, adjusting the time-axis of each 
of the template light curves by a common 'stretch' (s) factor (where 
the input SED time-series template is defined to have s = 1). The 
normalisations in the observed filters are independent, i.e., the ab- 
solute colours of the SED template used in the fit are not important 
and do not influence the fit. The result is a set of maximum-light, 
observer frame magnitudes, and a stretch factor. 

Once this observer-frame SiFTO fit is complete, it can be used 
to estimate rest-frame magnitudes in any given set of filters, pro- 
vided there is equivalent observer-frame filter coverage. This is 
performed by adjusting (or 'warping') the template SED to have 
the correct observed colours from the SIFTO fit, correcting for ex- 
tinction along the line of sight in the Milky Way, de-redshifting, 
and integrating the resultant rest-frame SED through the required 
filters. This process is essentially a cross-filter fc-correction, with 
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Rest-frame day 
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Figure 1. The light curve of SN2010jn/PTF10ygu in apparent AB magnitudes. The photometry has been taken from the Palomar 48-in telescope (P48) in R 
(solid red circles), the Liverpool Telescope (LT; solid green and orange circles) and Faulkes Telescope North (FTN; open green and orange circles) in g and 
r. No correction for extinction in the Milky Way or the SN host galaxy has been applied. Over-plotted is a SiFTO light curve fit (solid lines for P48 and LT 
data, and dashed lines for FTN data). Only data points up to +15 days past the (rest-frame) B-band maximum are used in the fit. The template lines are dotted 
after this phase and show the predicted evolution at later phase; however the SiFTO template in R/r is not expected to reproduce the observed secondary 
maximum. The arrows mark the epochs (given relative to the rest-frame B-band maximum) of the four HST spectra. 



the advantage that all the observed data can contribute to the over- 
all SED shape. We only use P48 and LT data (which have a wider 
phase coverage than the FTN data) for this colour warping; SiFTO 
does not handle filters that are very close in wavelength in this pro- 
cess (e.g., LT g and FTN g). However, all the FTN data are used to 
constrain the stretch and time of maximum light. 

Our SiFTO light curve fit is shown in Fig. Q] The fit is rea- 
sonable: the x 2 °f the fit is 27 for 43 degrees of freedom. Maxi- 
mum light in the rest-frame B-band was on MJD 55495.8 ± 0.1, 
or UT 2010-10-26.8, i.e., the SN was discovered some 15 days 
before maximum light, as expected from the analysis of the first 
Gemini spectrum. The stretch is 1.06 ± 0.02, corresponding to 
Amis (B) = 0.9 mag - the evolution of the light curve is slightly 
slower than for a typical SN la. We measure peak rest-frame abso- 
lute magnitudes at the time of B-band maximum light of Mb = 
-17.94 ± 0.04 and My = -18.36 ± 0.02, and a SiFTO colour 
C (broadly equivalent to B — V at B-band ma ximum l i ght) o f 
0.42 ± 0.03. We take the BV filter responses from lBesselll ll99Cl) , 
and these magnitudes are in the Vega system converted from the 
(near)-AB system in which the photometry is measured. The uncer- 
tainties are based on statistical errors propagated through the light 
curve fit. This indicates that PTFlOygu was fainter than its light 
curve stretch would indicate, and redder than a normal SN la, both 
presumably the result of extinction by dust (see Sec. |3.4| l. 

We can estimate the (B-band) rise time of the SN (r r , de- 
fined here as the time from explosion to B maximum) from the 
very early P48 data, independently of the SiFTO light curve fit. 
To this end, we measure the explosion epoch by fitting the P48 
R data at early epochs (< 5 d after discovery) using the analyti- 
cal equation f(t) = a x (r + T r ) 2 ('fireball model' ; iRiess et all 



1 19991 ; iGoldhaber et al.|[l998h . Here, r = t Q /(l + z) where t D is 
the (observer-frame) epoch relative to B maximum, and a is a nor- 
malising constant. This yields t t .r k = 18.6 ± 0.3 d. Leaving the 
exponent as a free fit parameter [f(t) — a x (r + r r )''], we obtain 
a best-fit rise time of r r ,f 1C c = 19.1 ± 1.2 d with an exponent of 
r\ = 2.3 ± 0.6 , slightly deviating from the value of rj = 1.8 derived 
bv |Pirol ( f20T2h . Enforcing his rj = 1.8 gives r = 18.3 ± 0.3 d. Our 
light curve template (stretch 1.0) has a rise time of ~ 17 d. Scaling 
it to the stretch of SN2010jn (i.e. to the higher opacity in such a 
slowly-declining object), we obtain an estimate of ryxcfit ~ 18.2d, 
slightly shorter than the values inferred from the early rise (and in- 
dicating the very early-time SiFTO templates may not match the 
data well). Our spectral models, on the other hand, lead us to prefer 
a somewhat higher value (~ 20 d; Sec. 13.41 Appendix [A), indicat- 
ing some tension between the rise time inferred from the early-time 
data with 77 = 2, and the rise time favoured by spectral modelling. 

3 SPECTRA AND MODELS 
3.1 Objectives 

With the aim of obtaining a comprehensive picture of SN2010jn, 
we now derive ejecta models from th e observed spectra. We use the 
method of Abundance Tomography JStehle et alj|2005h : an ejecta 
model is set up, where the density distribution is that of an ex- 
plosion model, but the abundances are left as 'free' parameters. 
Synthetic spectra are then calculated from this abundance / density 
model and the abundances are optimised so as to match the ob- 
servations with the synthetic spectra as closely as possible (cf. 
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Figure 2. W7 an d WS15DD3 (short WDD3) density profiles 
hwamoto et al.ll999T) used for our tomography experiments. 



Sec. 13.2.2b . Our code (Sec. 13.2. It was successfully applied to 
obtain spectral models for single UV spectra of various SNela 
dSauer et al.ll2008l : IWalker et alj|2012h . We profit from this expe- 
rience here, and perform for the first time tomography of a SN la 
including UV spectra, deriving a detailed abundance stratification. 
The UV spectra, especially those taken as early as — 10.5 or —5.8 d, 
when the outer layers leave their imprint in the spectrum, allow us 
to constrain iron-group abundances in the outer and intermediate 
layers of the SN. We perform the tomography for different explo- 
sion models. Thus, the quality of the spectral fits can be compared. 
Furthermore, for each explosion model we can check the consis- 
tency of the explosive nucleosynthesis, as calculated in the hy- 
drodynamic model, with the abundances we derive from the spec- 
tral fits. Finally, we can judge which hydrodynamic model offers a 
more realistic representation of SN 2010jn. 

Two different explosion models are tested in this pa- 



llwamoto et al. 


19991) and 


WS15DD3 i 


Iwamoto ct al. 



the 



'delayed-detonation' 
[). Both represent 
explosions. Their 



model 
single- 
density 



degenerate Chandrasekhar-mass 
distributions are shown in Fig.[2] The W7 model has been shown to 
provide excellent fits to light curves and spectr a of average as well 
as lower-than-avera g e luminosity SNIa (e.g. I Stehle et aI1l2005l : 
iMazzali et all 120081 : iTanaka et all 1201 ll) . Although its physical 
foundations are now known not to be realistic, as it assumes 
that a deflagration flame with parametrised propagation speed 
burns the star, its density distribution is representative of most 
explosion models, except for the outermost layers, so we use it 
as a benchmark. WS15DD3, a delayed-detonation model, is more 
physically consistent, and it has a higher explosion kinetic energy 
and 56 Ni mass, both of which turn out to be more appropriate for 
the (intrinsically) luminous SN2010jn. 

3.2 Methods 

We briefly describe the code which we use to calculate the synthetic 
spectra and the method of Abundance Tomography. 



3.2.1 Monte-Carlo radiative transfer code 



We u s e a Monte-Carlo radia ti ve transfer c ode I Abbott & Lucvl 
19851 IMazzali & Lucvl 1 19931 iLucvl 1 19991 IMazzalil |2000| and 
Stehle et al.l2005l) in spherical symmetry to compute the formation 



of the spectrum in an expanding SN envelope at a given epoch. 
The code assumes a sharp photosphere from which radiation is 
emitted, and computes the interaction of the photons with the ex- 
panding SN ejecta ('atmosphere'). An initial density profile, usu- 
ally the outcome of a hydrodynamic calculation, is used to cal- 
culate the densities within the atmosphere at the epoch required. 
This takes advantage of the fact that the ejecta are in a state of 
homologous, force-free exp ansion at the epochs relevant here (e.g. 
iRopke & Hillebrandj200a) . i.e., for each particle within the ejecta, 
r = vxt is a good approximation, where r is the distance from the 
centre, t the time from explosion and v the velocity of homologous 
expansion. Radius and velocity can therefore be used interchange- 
ably as spatial coordinates. 

The outward-flowing radiation at the photosphere, which is 
characterised by a velocity u p h, is assumed to be a Planck con- 
tinuum at a temperature T P h [/+ = B V (T V ^)}. Althou gh this 
is a relatively rough approximation teongardet"al]|2008h . espe- 
cially at late times an d in SNe where line opacity dominates (e.g. 
lHachinger et al.l20 09. and it can result in a mismatch in the overall 
flux in the red and infrared, it has only minor consequences on the 
determination of the abundances from prominent spectral lines. 

We simulate the propagation of the emitted photons through 
the envelope, considering 'photon packets' which undergo Thom- 
son scattering and bound-bound processes. The latter processes are 
treated in the Sobolev approximation, and a good approximation 
to the bound-bound emissivity is ensured by a downward branch- 
ing scheme. The constru ction of the simulation enforces radiative 
equilibrium l lLucyll 19991) . When packets are scattered back into the 
photosphere, they are considered to be re-absorbed and do not con- 
tribute to the output luminosity. An essential feature of the code is 
the consistent determination of a stationary state of the radiation 
field and of the excitation / ionisation struc ture of the plasma. To 
this end, a mod ified nebular approximation jMazzali& Lucvl 19931 : 
IMazzalil l2000h is used to calculate the gas state from a radiation 
temperature Tr and a dilution factor W cell per cell. The gas state 
and the radiation field are iterated in turn until the Tr values within 
the atmosphere are converged to the per cent level. The temperature 
at the photosphere, T p h, is automatically adjusted in these cycles so 
as to match a given bolometric output luminosity Lboi, compensat- 
ing for the re-absorption of packets. Finally, the emerging spectrum 
is not calculated from packet counts, which would introduce addi- 
tional 'Monte-Carlo noise', but rather via a formal int egral, using a 
source function from the packet statistics dLucvlll999l) . 



3.2.2 Abundance Tomography 

For an assumed density profile we model a series of spectra, infer- 
ring an optimum a bundance stratifica tion. This process of 'abun- 
dance tomography' l lStehle et al.l2005h is based on the idea that the 
opaque core of the expanding ejecta shrinks with time (in v space). 
In other words, the photosphere recedes, and deeper and deeper 
layers become visible and influence the spectral lines. 

From the earliest spectrum available we can obtain a photo- 
spheric velocity and abundances within the outer ejecta such that 
the synthetic spectrum matches the observed one. The next ob- 
served spectrum will still show absorption by the material in the 
outer zone, but additional layers, which only now are above the 
photosphere, also contribute to the absorption. Keeping the previ- 
ously inferred abundance values in the outer layers, the abundances 
in the newly exposed layers and a new photospheric velocity can be 
inferred. We continue this with later spectra. In some cases, abun- 
dances in the outer layers need to be revised in order to optimize 
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Figure 3. Spectra of SN2010jn (this work, black line) and SN1999ee 
iHamuY et al]|2002t orange / grey line) around maximum light. The spec- 
trum of SN 1999ee is displayed as observed, but has been multiplied by 
a factor to account for the difference in distance towards the two SNe 
(A/j, ~ 1.7 mag from NED - see acknowledgements). The spectrum of 
SN2010jn has been de-reddened by a total E(B - V) of 0.15 mag, and 
is then seen to coincide nicely with that of SN 1999ee. The spectral features 
discussed in the text are marked (features missing in SN 2010jn in grey with 
dashed marks). 



the fit to later spectra. In this case the entire spectral sequence is 
iterated in order to determine a consistent abundance distribution. 



3.3 Some general characteristics of SN 2010jn 

SN2010jn had a very broad light curve, corresponding to a high 
56 Ni mass, but it was a spectroscopi cally normal SNIa. It is 
in fact a close analogue of SN 1999ee lHamuv et alj J2002t> ; sec 
maxi mum-light optical spec tra in Fig. [3] rat her than of SN 199 IT 
(e.g. iFilippenko etai]|l992t) or SN 1999aa dGaravini et al]|2004l) : 
all of which, however, share similar light-curve properties. 

The UV flux in SN 2010jn is low, as in other SNe la, indicating 
significant line blocking. The optical P-Cygni absorption features 
(Fig. of Ca II (IR triplet) and Si II (A6355) show strong high- 
velocity feat ures (HVFs). Althou gh HVFs are common in SNe la at 
early times dMazzali etal]|2005h . those in SN2010jn stand out by 
their strength and high velocities (Tabled, similar to those found 
in SN 1984A dBarbon et alii 19891) . Only around B maximum do 
the HVFs vanish. 

Both Si III A4563 and Si II A5972 are very weak or absent 
in SN2010jn (Fig.[3]l. The Si III A4563 feature is characteristic for 
SNe la with high tem perature, while Si II A5972 is strongest in low - 
temperature objects jNugentetai]|l995l : lHachinger et ail 120081) . 
The weakness of both indicates a small abundance of intermediate- 
mass elernraits_(IME 2 _9_^ Z ^ 20), which is typical of luminous 
SNe dMazzali et al.l2007l) . Only the intrinsically strong Si II A6355 
line remains prominent. 



3.4 The rise time and reddening of SN 2010jn 

Since our code requires an epoch and a luminosity for each model 
spectrum as input, we need to make assumptions about the B-band 
rise time of the SN and the reddening. 

The B-band rise time r r (cf . Sec. 12. 21 is necessary to calcu- 
late the time from explosion t for each of the observed spectra. 
The spectral models presented here assume a rise time of 20 d. We 



have created tomography models for different values of the rise 
time (Appendix lAt . and find that a shorter rise time would im- 
ply higher (and probably unrealistic) iron-group abundances in the 
outer ejecta. This longer rise time implies that the SN after explo- 
sion has a 'dark phase' longer than predicted by the r 2 model. 

The reddening to SN 2010jn can be estimated in various ways. 
We have been able to measure the equivalent width of the host- 
galaxy Na I D line in our observed low-resolution spectra (2. 1 A on 
average). Thus, we can estimate the host-galaxy red dening using 
the m ore conservative (lower- reddening) relation of iTuratto et al .1 
d2003h : 

E(B - V) [mag] = EW(Na I D) [A] x 0.16 -0.01. 

Including galactic reddening, we arrive at a total reddening 
E(B - V^totai = 0.35 ma£| with this method. Since £W(Na I D) 
in low-res olution spectra may be a relatively uncerta in proxy for 
extinction dBlondin et al.l2009c|Poznanski et al.l2011D we obtained 
further independent estimates. 

We can infer the reddening of SN2010jn relative to a bright 
SN with similar lig ht curve shape, SN 1999ee (Ami 5 (B)= 0.95; 
lHamuv et alj 120021) . Comparison of the spectrum of SN2010jn 
at — 0.2d and a maximum-lig ht spectrum of SN 1999ee (Fig. [3j 
suggests a red dening difference of AE(B — V) tot ni ~ 0.15 mag. 
ISasdellil d201ll) modelled SN 1999ee in detail and found that it was 
most probably reddened by E(B — V)tom = 0.26 mag. Assum- 
ing that the SNe are almost identical, this gives a reddening of 
E(B - VOtotai ~ 0.41 mag for SN2010jn. 

Finally, our light curve fits give a peak optical colour of B — 
V — 0.42 mag (Sec. |2.2t , which can also be taken as an indicator 
of extinction. 

As the three extinction estimates are reasonably consistent 
with one other, we use their average, 

E(B-V) ta u = 0.39 mag. 



3.5 Models based on the W7 density distribution 

We begin by computing synthetic spectra using W7, in order to 
compare the resul ts for SN20 1 0jn wi th those obtai ned for other 
SNe [SN 2002bo, IStehle et alj d2005h: S N2003du. iTanaka et alj 
d201 lh : SN 2004eo7 lMazzali et alj d2008l) l. Since SN2010jn was 
more luminous than all these SNe, differences in the results may 
be expected. 

Our best-fit synthetic spectra for SN2010jn based on the W7 
model are shown in Fig.|4]together with the observations. We first 
discuss the spectral fits, then we point out some major shortcomings 
of the W7-based model. Finally, we show and discuss the abun- 
dance stratification we infer using the W7 density. 



3.5.1 Spectral models 

The first spectrum we model has an epoch of 12.9 d before B-band 
maximum, i.e. 7. 1 d after explosion (for t r — 20 d). The model lu- 
minosity is Lboi = 0.7 x 10 9 Lq. The spectrum shows all lines char- 
acterising a typical SN la in the photospheric phase. This indicates 
that not only O, but also IME (~33% by mass) and some iron are 



4 For low redshifts, such as that of SN 2010jn, galactic and host extinction 
values can simply be added up to obtain the total reddening. 
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Figure 4. l()jn-W7 model sequence, based on the W7 density profile (red lines). The observed spectra are shown in black for comparison; insets show the UV 
in more detail. Observations and models are shown de-redshifted (i.e. in the SN's rest frame); no correction for host-galaxy reddening has been applied to the 
observed spectra [i.e. the models have been reddened by E(B — V)h os t = 0-36 mag]. For the most prominent UV features, the respective ions dominating the 
absorption in our model are given (in order of importance). Note that Cr and Ti cause significant quasi-continuous absorption by a large number of weak lines, 
which are not reported here. 
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Figure 5. Abundances of W7 nucleosynthesis calculations jlwamoto et al.ll99^ top panel, plotted in mass space) vs. Abundance Tomography of SN 2010jn, 
based on the W7 density profile (10jn-W7 model - middle panel: in mass space; lower panel: in velocity space). The C abundance in our model is an upper 
limit, as no C feature is present in the observations. The Ni/Co/Fe abundances are given in terms of the mass fractions of 56 Ni and stable Fe at t = 
[X ( 56 Nio), X(Ferj)]; in our spectral models, no stable Ni or Co and no radioactive Fe are assumed to be present. 
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Table 2. Velocities of HVFs in SN2010jn/PTF10ygu. The values given for each feature correspond to the blueshift of the point of deepest absorption. While 
in S II A5640, this point normally corresponds to absorption near the photosphere, in Si II A6355 and the Ca II IR feature it is determined by a high-velocity 
component before maximum. In the later spectra, a photospheric and a HVF component can be distinguished in Si II A6355. 



Calendar date 


rest-frame phase 


v(Ca II A ~ 8600) 


v(Si II A6355) 


v(S II A5640) 


(UT) 


(days) 


flans- 1 ) 


(kms- 1 ) 


(km s- 1 ) 


Oct 13 


-12.9 


35000 


23000 


19000° 


Oct 16 


-10.5 


28000 


20000 


18000° 


Oct 20/21 


-5.8 


22000 


18000 


13000° 


Oct 26 / 27 


-0.3 


19000 


17500/14000° 


12000° 


Oct 31 / Nov 1 


+4.8 


17000 


17000/13500° 


11500° 



"values refer to photospheric absorption 



present above the photosphere, which is located at 14600 km s _1 
in our model. The observed lines show high blueshifts, which can- 
not be reproduced by the W7-based model as it lacks high- velocity 
material. The absorption troughs at ~ 4300A and ~ 4700A, which 
are both due to a mix of Fe and IME lines, serve as diagnostics for 
the Fe content in the absence of UV dat^f]. 

At 10.5 d before B-band maximum, the W7 -based model has 
a bolometric luminosity of Lboi = 1.8 x 1O 9 L0 and a photo- 
spheric velocity of 13350 km s~ . The temperature of the photo- 
spheric black-body spectrum has risen from T p h = 11160K (at 
- 12.9 d) to T ph = 12870 K. The low UV flux is indicative of effi- 
cient line blocking, and requires significant iron-group abundances 
in the photospheric layer [X( 56 Ni ) ~ 4%, X(Fe ) ~ 8% and 
X(Ct) ~ 0.1% at the photospherefl Reverse-fluorescence pro- 
cesses, induced by iron-group elements and I ME in the outer lay- 
ers, actually enhance the UV flux in SNela IjVlazzali 20o3)- m 
higher-luminosity models such as the ones for SN2010jn, how- 
ever, Fe-group elements also have a str ong absorbing effec t in the 
UV except for some opacity windows jWalker et al. |2012|) . Addi- 
tionally, the radiation field in these models is bluer, disfavouring 
reverse fluorescence. The two prominent UV absorption features in 
the - 10.5 d model (~3000 A and ~3200 A) are mostly due to Fe II 
and Ni Il/Co II lines, respectively. At < 2900 A, the UV flux is de- 
pressed by Fe II and Mg II lines (and other lines further bluewards). 
As the ~ 3200 A feature is dominated by Ni II and Co II lines, 
we could separate the effects of Fe and Ni / Co on the spectrum. 
Also Cr and Ti have distinct effects: as long as their abundances are 
reasonably low, they suppress the UV flux without making strong 
features. Thus it was possible to determine the abundance of each 
element separately. 

Five days later (— 5.8d, i.e. ~14.2d after explosion), the SN 
had brightened significantly (Lboi = 3.6 x 10 9 Lq), and the spec- 
trum had become somewhat bluer (as indicated by the higher T p h = 
13170 K). The photosphere in our model recedes to 1 1400 km s -1 , 
and is located in a zone which is dominated by Si, but where 
also iron-group material is abundant. The UV feature at ~3000 A 
now has a stronger contribution of Fe III (and a large contribu- 
tion by Si III), as the higher luminosity increases ionisation. Ad- 
ditionally, the feature at ~3200 A, which was dominated by Ni II 
and Co II lines in the earlier spectrum, is now contaminated by 



5 The ~ 4700A trough is slightly too deep in the synthetic spectrum at this 
epoch; however, it is slightly too shallow at the next epoch. 

6 The Ni / Co / Fe abundances are given in terms of the mass fractions of 
56 Ni and stable Fe at t = [X( 56 Ni ), X(Fe )]; no stable Ni or Co and 
no radioactive Fe are assumed to be present, as these species - assuming 
realistic abundances - have less of an impact on our models. 



Fe III. Also, strong Co III lines begin to contribute there. Because 
of strong back-warming effects, Si is highly ionised and little Si II 
remains. Therefore, the Si I I A5972 absorption (cf. iNugent et al.l 
ll995l : lHachinger et alj|2008h is weak. The O I A7773 feature is es- 
sentially missing for the same reason. We can also set rather strin- 
gent limits on the C content in the outer ejecta, as even small C 
abundances produce a visible C II A6580 line (which is not seen in 
the observations). 

For epochs at B-band maximum and later, SN2010jn dis- 
played a normal SNIa spectrum. Our models have a bolometric 
luminosity L bo \ = 4.6 x 10 9 L Q at -0.3 d and L bo i = 4.0 x 1O 9 L 
at +4.8 d, respectively. The HVFs in Si II and Ca II have become 
much weaker, as it is the case in most SNela. The model there- 
fore fits the observations better and better with time. Even though 
the HVFs disappear, line velocities remain high compared to other 
SNela. The photospheric velocity is ~ 8850km s _1 even around 
B maximum. At +4.8 d, the photosphere is at ~ 7200 km s -1 , 
and the photospheric temperature has dropped from from T p h = 
13450 K (the value at -0.3d) to T ph = 1251 OK. High iron-group 
mass fractions (> 85%) are needed at 7200 < v < 11400 km s _1 
in order to fit the low UV flux after maximum. These velocities 
are high, indicating efficient 56 Ni production and a luminous SN. 
The UV features after maximum are dominated by lines of Fe II 
(~ 3000 A) and Co Ill/Co II (~ 3200 A). Flux blocking around 
~2900 A still is mostly due to Fe II and Mg II lines. 



3.5.2 Shortcomings of the W? '-based model 

The major problem of the W7 -based model is that the lines in the 
synthetic spectra at the earlier epochs generally do not reach ve- 
locities (blueshifts) as high as the observed features, although we 
adopted photospheric velocities as high as possible. This is a result 
of the steep decline of the W7 density at high velocities. The mis- 
match can be improved by assuming a delayed-detonation model 
(Sec. 13. 6t . We did not attempt to fit the HVFs, at this would re- 
quire the density and /or electron density in the outermost zone 
to be enhanced prob ably by mixing with a circumstellar medium 
dAltavilla etai1l2007l) . which is beyond the scope of this paper. 

Another notable shortcoming of the models is a lack of flux 
around 4000A, in particular at — 10.5 d. The high flux in the ob- 
served — 10.5 d spectrum may not be completely reliable (there is a 
sudden change with respect to — 12.9 d), but more absorption in the 
Ca H&K feature would increase the re-emission into the peak. 

Later synthetic spectra tend to show an excess of flux in the IR 
(redwards of ~6500 A) compared to the observations. This is rel- 
atively independent of the assumed density model (cf. the WDD3- 
based models in Sec. |3.6| l. It is rather an artefact of using a black- 
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body spectrum at the photosphere (Sec. 13.2.7V The resulting mis- 
match is however not expected to have major effects on our abun- 
dance determination. The fact that this happens so early suggests 
that the black body approximation fails earlier than in most SNe la. 
This is probably the consequence of the presence of significant 
amounts of 56 Ni near and above the photosphere even when the 
photospheric velocities are still high. 



3.5.3 Abundance structure 

The spectra probe the ejecta structure from the outermost layers 
down to 7200 km s _1 (+4.8 d photosphere). The abundances are 
shown and compared to the W7 nucleosynthesis in Fig. [5] 

The model contains two zones with reduced abundances of 
burning products at 22000 km s _1 ^ v ^ 33000 km s _1 and v > 
33000 km s _1 . These are not strictly needed to improve the fit to 
the spectra, but were added in order to have a structure analogous 
with the WDD3-based models (see Sec. |3.6l l. 

Moving inwards, there is a zone extending over ~ 0.2 Mq 
(14500 km s" 1 < v < 22000 km s _1 ) which is dominated by O, 
but already contains a significant fraction of IME (Mg, Si and S) 
and iron-group elements. 

IME clearly dominate in the velocity range 11000 km s" 1 < 
v < 14500 km s _1 (i.e. over a mass range of ~ 0.2Mg). Below 
~ 11000 km s _1 (at a mass coord inate of > 0.9Mq ), where in 
most SNe la IME are still dominant l lMazzali et alj2007l) . the ejecta 
of SN2010jn begin to be dominated by iro n-group material. Th is 
has only been seen in very luminous SNe la jMazzali et al.ll2007b . 

Comparing our tomography result with the W7 nucleosynthe- 
sis calculation (Fig.O reveals that iron-group elements are present 
in SN2010jn at higher velocities. The higher burning efficien cy is 
consistent with the luminosity (see also IStritzinger et alj2006l) and 
also with a high kinetic energy, as indicated by the high velocities 
in all spectral lines. Also wit h respect to other SNe studied with the 
tomography technique (e.g. IStehle et al. II 20051 ; iTanaka et al .11201 ll ; 
lMazzali et 'ai]|200a iHachinger et alj|2009l) . the SN shows a far 
more efficient burning. An additional difference is that SN 2010jn, 
unlike these other SNe, exhibits rather sharp transitions between the 
zones dominated by different elements (Fe-group, IME, C/O). This 
may have to do with the fact that in an efficiently burning SN la the 
zone in which IME are produced and the zone where O is left are 
constrained to relatively narrow outer shells, and less "mixing" is 
possible. 



3.6 Models based on the WDD3 density 

Having noticed the shortcomings of W7 as a backdrop density pro- 
file, we now test whether the spectra of SN2010jn can be better 
fitted using the density of a delayed-detonation model. In a delayed 
detonation, an initial subsonic burning phase is followed by a phase 
of supersonic burning which more effi ciently converts the original 
C/O mixture into processed elements l lKhokhlovIl 1 99 ll) . This type 
of explosion can produce a larger amount of iron-group elements 
as well as IME, and a higher kinetic energy, placing more mate- 
rial at high velocities. Although this scenario appears to be very 
promising for SNe la, the details of the transition from subsonic to 
supersonic burning are unfortunately still unclear. 

Given the high luminosity of SN 2010jn and the efficient 
production of Fe group elements we find in our model, we use 
model WS15DD3 (in short WDD3), the brightest and most en - 
ergetic of the WS15DD model series of llwamoto et alj Jl999h . 



Table 3. Luminosities, photospheric velocities and temperatures of the pho- 
tospheric blackbody in the 10jn-WDD3 model. 



Rest-frame phase 


Lbo\ 


I'ph 




(days) 


(1O 9 L ) 


(km s" 1 ) 


(K) 


-12.9 


0.75 


14600 


10220 


-10.5 


1.7 


13350 


11670 


-5.8 


3.5 


11400 


12830 


-0.3 


4.5 


8850 


13430 


+4.8 


4.0 


7200 


12700 



WDD3 synthesizes 0.77 Mq of Ni and has an explosion energy 
E km = 1.43 foe. 

For spectral modelling, the most important difference in the 
density structure of WDD3 in comparison to W7 (cf. Fig.|2j is the 
higher density in the outer part of WDD3. 

In Fig. 1 our WDD3-based optimum models for SN2010jn 
are compared to the data and to the W7-based models presented in 
the previous section. 



3.6.1 Spectral models 

The synthetic spectra we obtained using the WDD3 model repro- 
duce the high line velocities of SN2010jn significantly better [ex- 
cept for the HVF, wh i ch we do not attempt to reproduce - cf. Sec. 
l331and lTanaka et alj j20lTI) l. At the earliest epochs, the larger ex- 
tent of the line-forming zone in WDD3 (towards high velocities) 
makes the WDD3-based spectra smoother, which matches the ob- 
servations much better. In particular, for the two earliest epochs 
(-12.9d, -10.5d), the fit to Si II A6355 and the Fe-dominated trough 
at ~4700A are broader and give an improved fit over W7. The 
Ca II H&K features shows a significant improvement even at later 
epochs, and the flux at 4000A matches the observations somewhat 
better. The UV region is matched just as well as in the W7-based 
model. 

The improvements with respect to W7 are basically due to 
the higher density in the outer zones in the delayed-detonation 
model. The increased amount of material at high velocities en- 
hances the absorption at larger blueshifts directly, but also through 
an indirect effect: the higher densities lead to a decreased ionisa- 
tion, favouring singly-ionised species over doubly-ionised ones. As 
many important lines in SNIa spectra are due to singly-ionised ma- 
terial, the blue line wings are better reproduced. However, apart 
from the outer layers, the changes in the model atmospheres from 
W7 to WDD3 are limited. This reflects in the fact that the op- 
timum WDD3 model could be constructed using practically the 
same values for photospheric velocity and bolometric luminos- 
ity as for W7 (cf. Table [3). The photospheric black-body tem- 
peratures are similar to the W7-based models except for the first 
two epochs, where the WDD3-based models have a flatter tem- 
perature gradient because WDD3 has a lower density than W7 at 
12500 km s" 1 < v sC 16000 km s" 1 (cf. Fig.|2j. 

The UV features in the WDD3-based model at the earliest 
epochs are dominated by the same elements as in the W7-based 
model (Fe near ~ 3000 A, Ni, and Co near ~ 3200 A. At later 
times, all Fe-group elements contribute somewhat more to all fea- 
tures because of increased line opacities, but still the effects of the 
different Fe-group species are sufficiently distinct that the abun- 
dance of each element can be constrained. 



© 2012 The Authors. Journal compilation © 2012 RAS, MNRAS 



Analysis of SN 2010jn ultraviolet spectra 1 1 



c 

ZJ 

1 

CO 
1 

i 

CO 





| i i i | i i i | i i i | i i i | i i i | i i i | i i i | i 

WDD3-based model, — SN 201 Ojn: -12.9d 



t=7.16d, lg(L/L,. : )=8.88, v Dh =14600km/s 



I 1 1 1 I 1 1 1 I 1 1 1 I 1 1 1 I 1 1 1 I 1 1 1 I 1 1 1 I 1 1 1 I 1 1 1 

WDD3-based model, — SN 201 Ojn: -10. 5d 

t=9.46d, lg(L/L,,)=9.23, v ph =13350km/s 

W7-based model 




3500 X[A] 



WDD3-based model, — 
t=14.20d, lg(L/L,,)=9.54, v ph =1 1400km/s 

W7-based model 



SN 201 Ojn: -5.8d 




WDD3-based model, — 
t=19.66d, lg(L/L @ )=9.65, v ph =8850km/s 

W7-based model 




SN 201 Ojn: -0.3d 







L- ' ' \ 






O 


o 






o 


O 




CD CD 


CD = 






LL LL 


LLC/5 




i 











3000 



+ 



3500 A.[A] 



■ I ■ ■ ■ I ■ 



WDD3-based model, — SN 201 Ojn: +4.8d 

|l=24.78d, lg(L/L l£; )=9.61 , v ph =7200km/s 

W7-based model 






■_»■ ■ i ■ ■ 


. °> o 


o 


S o 


O 


CD CD CD 




LL LL LL 

: — t^WT 





3000 



3500 A.[A] 



3000 



5000 



7000 



9000 



11000 



MA] 



Figure 6. 10jn-WDD3 model sequence, based on the WDD3 density profile (light, blue lines). Observed spectra (black lines) and 10jn-W7 spectra from Fig. [4] 
(red, dotted lines) are given for comparison; insets show the UV in more detail. The differences in the density profile influence the line velocities. Observations 
and models are shown de-redshifted (i.e. in the SN's rest frame); no correction for host-galaxy reddening has been applied to the observed spectra [i.e. the 
models have been reddened by E(B — V)h os t = 0.36 mag]. For the most prominent UV features, the respective ions dominating the absorption in our model 
are given (in order of importance). Note that Cr and Ti cause significant quasi-continuous absorption by a large number of weak lines, which are not reported 
here. 
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3.6.2 Abundance structure 

In Fig. [7] we show the abundances of our 10jn-WDD3 model to- 
gether with a plot of the original WDD3 nucleosynthesis. This 
model may be regarded as a reference for SN2010jn: the good 
fits to the observed data indicate that the density of WDD3 pro- 
vides a reasonable description to the actual ejecta. The abundance 
stratification of the WDD3-based models in velocity space is qual- 
itatively similar to that of the W7-based models d3.5.3t . despite the 
differences in the explosion models. Therefore, we can say that our 
results on the chemical structure should be reasonably independent 
of the density model. 

In order to construct our optimum models, we have intro- 
duced two independent abundance zones in the outermost part of 
the ejecta (22000 kms'^K 33000 km s _1 and v ^ 33000 km 
s _1 ). These zones contain fewer burning products, ensuring op- 
timum spectral fits (e.g. by avoiding spurious high-velocity ab- 
sorption in the Fe/Mg feature at ~ 4300 A). The Ca, Si and S 
abundances at 22000 km s _1 ^ v ^ 33000 km s _1 have been opti- 
mised to fit the higher- velocity parts of the respective features (ex- 
cept for the HVFs). Despite the decrease of most burning products 
(with respect to the lower layers), Ca and Ti / Cr remain abundant, 
possibly the result of incomplete burning. 

Our model yields an upper limit of ~60% to the combined 
abundance of C and O at 14600 km s _1 < v < 33000 km s^ 1 . 
With higher mass fractions, we would have lower fractions of burn- 
ing products, and this would degrade the fit. 

The moderate C/O abundance between 14600 and 33000 km 
s _ is compensated by an IME mass fraction of >30%, a mass frac- 
tion of directly synthesised Fe of ~3%, and a mass fraction of ~2% 
of 56 Ni and decay products. The UV spectra make it possible to 
disentangle the effects of Fe, 56 Ni and Ti/Cr, and to infer their 
respective abundances (see Fig. [7}. 

Assuming the WDD3 density profile, we have surveyed 
0.8M Q of the ejecta (down to the +4.8 d photosphere at 7200 km 
s _1 ). The C/O content in this region is ~ 0. 15 Mq. The mass 
in IME is ~0.21 Mq, and ~O.45M0 are in iron-group elements 
(~0.42M Q in 56 Ni). 

If the ~ O.6M0 core, which cannot be analysed with our 
spectra, has a 56 Ni mass fraction as at our innermost photosphere 
(90%), SN 2010jn contains 0.9 M of 56 Ni. This is consistent with 
a very luminous SN la. In order to complete our anal ysis at the low- 
est ve locities, nebular spectra would be required (cf. lMazzali et al.l 
l20Tlh . 

Finally, comparing our abundance diagram with the delayed- 
detonation nucleosynthesis (Fig.UJ, we find good agreement in the 
general structure with few exceptions. SN2010jn shows somewhat 
more mixing than WDD3: in the IME-dominated zone, both oxy- 
gen and (to a smaller degree) iron-group elements are abundant. 

A main finding of this analysis is the presence of burned ma- 
terial in the outer layers of SN2010jn. The high iron-group abun- 
dances (Fe, Ti/Cr, 56 Ni) in these zones, which lead to a high total 
opacity of the ejecta, explain the large stretch of SN2010jn. Evi- 
dence for 56 Ni out -mixing has already been given for the less en- 
ergetic SN2011fe iNugent et al1l201lMPirdl2012h , using different 
diagnostics. 



3.7 The effect of iron-group material in the outer layers on 
the spectra 

The determination of the abundances of different iron-group ele- 
ments in the outer layers of SN2010jn (Sections l3.5N3.6t depends 
crucially on the early-time UV spectra. 

To demonstrate this, we computed spectra using modified ver- 
sions of our WDD3-based ejecta model (Fig. [8j: we changed the 
iron-group abundances up and down by a factor of two in the outer 
layers above the — 10.5 d photosphere (13350 km s _1 ), where the 
ejecta are dominated by O and IME (cf. Fig.[7j- 

The earliest spectra (— 12.9 d, —10.5 d) show a clear sensitiv- 
ity of the UV flux on iron-group content, while changes in the op- 
tical are rather moderate. This is due to the fact that iron-group el- 
ements - apart from Fe itself - do not produce pronounced absorp- 
tions in the optical. Even the optical Fe lines are strongly blended 
with lines of other elements (Mg, Si). Therefore, the UV is the only 
reliable quantitative diagnostic of iron-group content in the outer 
ejecta. 

At — 5.8d, the UV still reacts, but somewhat less, since the 
photosphere is located deep inside the region affected by the 
changes. The optical Fe lines again do not seem to react to metal 
content. The change seen in the Si II A6355 region is in fact due to 
a change in back-warming, which is however only an indirect and 
crude diagnostic for metal abundance. 

At later epochs the UV spectra are more sensitive to the metal 
content in the inner layers. They remain practically unchanged 
when the outer zones are modified (Figure [8] lower two panels). 
They may be useful to diagnose the iron-group content in the core, 
which however can be better analysed using nebular-phase spectra. 



4 SUMMARY AND CONCLUSIONS 

We have used a well-established spectrum synthesis code to 
conduct, for the first time, Abundance Tomography with a se- 
ries of photospheric UV/ optical spectra of a SNIa, SN2010jn 
(PTFlOygu) The early observations have proved extremely valu- 
able in order to constrain the iron-group content in the outer lay- 
ers. Early UV spectra are probably the best diagnostics for this, as 
we have shown by recalculating our original models with modified 
iron-group abundances. The later spectra allow us to constrain the 
abundances further inwards and to obtain a comprehensive picture. 

From our abundance analysis, we find that SN2010jn was an 
energetic explosion with significant amounts of intermediate-mass 
and iron-group elements in the outer layers (v > 13350 km s _1 ). 
SN2010jn has only a thin oxygen-dominated zone and a limited 
IME zone. Even these layers contain directly synthesised Fe (~ 
0.01 M Q ) as well as 56 Ni and decay products (also ~ 0.01 M Q in 
total). The presence of iron-group elements leads to a high opacity 
of the SN, consistent with the slow light-curve evolution. The UV 
spectra allowed us to determine the content of iron-group elements 
and to differentiate among them. 

We favour for SN20 10jn a Chandrasekh ar-mass delayed- 
detonation model (WDD3. Ilwamoto et al.l fl999l) . It allows us to 
reproduce the high expansion velocities in the observed spectra 
better than the 'fast deflagration' model W7 jNomoto et ai1(l984l : 
Ilwamoto et al.fl999T ). which has little material at high velocities. 

The presence of Fe and 56 Ni in the outermost layers is not pre- 
dicted by the original (1-D) WDD3 model. The abundance of Fe in 
the outer layers (13350 km s _1 ^ v ^ 33000 km s _1 ) is between 
one and two orders of magnitude above solar. This can not come 
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from the progenitor; it is rather a clue about the explosion proper- 
ties of SN 2010jn. A high Fe-group abundance is consistent with an 
energetic SN, and the opacity due to Fe and 56 Ni in the outer lay- 
ers explains the small decline rate of SN2010jn. It is a challenge 
for explosion models to reproduce this. Multi-dimensional models 
may show some dumpiness or asymmetry which could appear to 
an observer as 'outwards-mixing' of iron-group material. 

Our models favour a longer rest-frame rise time (~20 d) than 
predicted by the empirical r 2 analysis of the early light curve 
(18.6 d). A longer rise time would imply that the SN after explo- 
sion has a 'dark phase' longer than predicted by the r 2 ('fireball') 
model. SN 2010jn is more luminous than average SNe la: it remains 
to be seen whether the same phenomenon is found in other objects. 
Very early data are best suited to tackle this issue. 

Finally, we remark that a precise analysis of the abundances in 
the outermost layers, as presented here, should be complemented 
by an analysis of the ejecta core, for which nebular spectra are nec- 
essary. Also, observations of more SNe la in the UV will be needed 
to establish the properties of the outermost layers of SNe la with 
different characteristics and light-curve properties. 
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Figure Al. Distribution of the most abundant / relevant elements (in veloc- 
ity space) in our T r = 18.6d models compared to the those in the standard 
(r r = 20 d) models. Top panel: W7-based model; bottom panel: WDD3- 
based model. The black and violet graphs [mass fractions of 56 Ni and sta- 
ble Fe at t = 0: A"( 56 Nirj), X(Feo)] give the abundances of the dominant 
iron-group elements, which are unrealistically high in the outer layers of the 
T r = 18.6d models (see text). 



APPENDIX A: CHOICE OF THE B-BAND RISE TIME 

In order to calculate the models presented in the main text we have 
adopted a value of 20 d for the B-band rise time (in the rest frame 
of the SN). Here we discuss why we have chosen this rise time, 
which is a bit longer than the one which is inferred from the ob- 
servations assuming a t 2 behaviour for the rising part of the light 
curve (r r ,flre = 18. 6d; Sec. 12.2b . 

During the process of modelling the SN, we have created vari- 
ous test models in order to determine which rise time allows for op- 
timum models in terms of fit quality and physical consistency. Here, 
we show the differences between models assuming T r ^ IC — 18.6 d 
and t t = 20 d, which lead us to accept the latter value as the one 
more compatible with our models. 

The fundamental difference when going from models with 
t t = 20 d to models with a shorter rise time (r r = 18.6d) is 
the smaller size of the ejecta. This affects the expansion state of 
the ejecta, but also the radius of a photosphere located at a certain 
velocity v p h- This affects the spectrum at the photosphere: a higher 
temperature T p h is required to emit the same amount of radiation 
from a smaller photosphere. 

To have photospheres at the same radius, keeping a similar 



colour and temperature as in the r r = 20 d model, a higher v p h 
must be used if r r = 18.6 d. For the earliest spectra (— 12. 9 d, 
— 10.5 d), Vph in our standard (20 d) model is already very high. 
An even higher velocity would mean that almost no mass is above 
the photosphere, but then line absorption would be reduced and it 
may then be difficult to reproduce the low observed UV flux (i.e. 
the flux in the model will be too high). The alternative would be to 
preserve the photospheric velocity from the standard calculation. 
In this case, however, the photospheric spectrum would be bluer 
(higher T p h required to emit the same luminosity from a small- 
radius photosphere, as mentioned above). Again, one would be left 
with too much UV flux. 

The only way to circumvent this and fit the observed spectra 
with t t = 18.6 d is to use higher mass fractions of the iron-group 
elements in the outer ejecta in order to block the UV flux. Our stan- 
dard (T r = 20 d) models have already shown a relatively strong 
out-mixing of iron-group material (Sections 13.5.31 13.6.2t with re- 
spect to the orig inal W7/WDD3 models dlwamoto et al J 19991) . As- 
suming T r — 18.6 d, the iron-group fractions in the outer layers are 
even higher and less realistic. The most relevant abundances in the 
t t — 18.6 d and r r = 20 d models are compared in Fig. I A 1 1 In 
the W7 -based model with r r = 18.6 d, an iron-group abundance 
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(sum of all iron-group elements) of ~15% is needed even above 
15000 km s _1 . This is about twice the result of the explosiv e nucle- 
osynthesis in the original W7 model Jlwamotoet all 19991) . which 
constitutes an inconsistency. In general, it will be difficult for most 
explosion models to produce high abundances of Fe-group mate- 
rial in the outermost layers: the densities there are small, leading 
to incomplete burning, and only "mixing" can bring burned ma- 
terial outwards. Our WDD3-based model shows a smaller effect, 
but still the high iron-group abundance in the outer layers of the 
T r — 18. 6d model (~ 7-23%, with a stronger gradient than in the 
W7 case) leads us to prefer a rise time of > 20 d, for which the 
iron-group abundances in the outer layers of our models seem con- 
sistent with the WDD3 results if one assumes moderate additional 
out-mixing (cf. Sec. 13.6.21 . 

Models for UV observations at even earlier epochs would re- 
act even more sensitively to the assumed rise time. If earlier UV 
observations were available for SN2010jn, a stricter lower limit on 
the rise time may have been determined through modelling. In or- 
der to exploit this possibility for future objects, we encourage UV 
observations of SNe la at epochs as early as possible. 
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